doseresNMA antidep/data.addFluxdose.R

# add fluoxetine-equivalent doses
data_with_fluxdose <- function(data=antidep){
  data$dose_ratio <- as.numeric(with(data, ifelse(
    drug=="agomelatine",0.75,ifelse(
      drug== "amitriptyline",0.33,ifelse(
        drug=="bupropion", 0.11, ifelse(
          drug=="citalopram",1,ifelse(
            drug=="clomipramine",0.34,ifelse(
              drug=="desvenlafaxine",0.40, ifelse(
                drug=="duloxetine",0.33, ifelse(
                  drug=="escitalopram",2.22,ifelse(
                    drug=="fluoxetine",1,ifelse(
                      drug=="fluvoxamine",0.28, ifelse(
                        drug=="levomilnacipran", 0.67,ifelse(
                          drug=="milnacipran",0.2, ifelse(
                            drug=="mirtazapine",0.79,ifelse(
                              drug=="nefazodone",0.07,ifelse(
                                drug=="paroxetine",1.18,ifelse(
                                  drug=="reboxetine",3.48,ifelse(
                                    drug=="sertraline",0.41,ifelse(
                                      drug=="trazodone",0.1,ifelse(
                                        drug=="venlafaxine",0.27,ifelse(
                                          drug=="vilazodone",2,ifelse(
                                            drug=="vortioxetine",2,"0"
                                          )
                                        )
                                      )
                                    )
                                  )
                                )
                              )
                            )
                          )
                        )
                      )

                    )
                  )
                )
              )
            )
          )
        )
      )
    )
  ))
  )
  data$dose <- data$dose_ratio*data$dose
  return(data)
}
#
htx-r/doseresNMA documentation built on Jan. 28, 2021, 5:32 a.m.